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ABSTRACT 



Aims. We simulate the light variability of the Ap star e UMa using the observed surface distributions of Fe, Cr, Ca, Mn, Mg, Sr, and 
Ti obtained with the help of the Doppler imaging technique. 

Methods. Using all photometric data available, we specified light variations of s UMa modulated by its rotation from far UV to IR. 
We employed the LLmodels stellar model atmosphere code to predict the light variability in different photometric systems. 
Results. The rotational period of e UMa is refined to 5^088631(18). It is shown that the observed light variability can be explained as a 
result of the redistribution of radiative flux from the UV spectral region to the visual caused by the inhomogeneous surface distribution 
of chemical elements. Among seven mapped elements, only Fe and Cr contribute significantly to the amplitude of the observed light 
variability. In general, we find very good agreement between theory and observations. We confirm the important role of Fe and Cr in 
determining the magnitude of the well-known depression around 5200 A by analyzing the peculiar a-parameter. Finally, we show that 
the abundance spots of considered elements cannot explain the observed variabilities in near UV and /? index, which probably have 
other causes. 

Conclusions. The inhomogeneous surface distribution of chemical elements can explain most of the observed light variability of the 
A-type CP star s UMa. 
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1. Introduction 

Chemically peculiar (CP) stars are main-sequence A- and B-type 
stars that have a number of characteristic photometric and spec- 
troscopic anomalies. Slow rotation (and, consequently, weak 
meridional currents), globally structured magnetic fields, and 
weak convection allow a build up of the prominent abundance 
anomalies and inhomogeneities as a result of p article diffusio n 
driven by the radiation field (first proposed by Mic haudlll970l) . 
These inhomogeneities are frequently s een as abundance spots 
on stellar surfaces (see, for example, iKhokhlova et al.l l2000t 
Kochukh ov et al.ll2004t iLehmann et al.ll2007l) . and as clouds of 
chemical ele ments concentrated at di f ferent heights in stellar 
atmospheres fe vabchik ova et al.ll2002l l2008t iKochukhov et all 
2006; Sh ulyak et al.ll2009l and references therein). 

Among the defining characteristics of CP stars, their light 
variability in different photometric systems still awaits quanti- 
tative explanation. That this variability is seen in independent 
systems, simultaneously in narrow- and broad-band bands, and 
in anti-phase in far-UV and optical spectral regions suggests 
that the global flux redistribution caused by the phase-dependent 
absorption is t he pr obable cause of the observed phenomena 
(Mol nanll973L I1975I). If so, there may be a point (or points or 
even regions, see Mikulasek et al. 2007b) in the spectra where 
the flux remains almost unchanged, and vary more or less in anti- 
phase on both sides from it. This region, called a "null wave- 
length region", has been noted and confirmed by a number of 



studies (see, for example. 



Leckrond[l974t iMolnaii [19731: Hamarl 
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1 1 9771: ISokolovll2000l 120061) . However, that the truly "null wave- 
length irejrijjn" inj^rs_with more complex light curves may not 
exist (Sokolovl l2TMl2TjIol) , The connection between the abun- 
dance spots and flux redistribution has emerged as a natural ex- 
planation of the rotationally modulated light curves of CP stars, 
but the details of this relation still have to be determined on a 
solid theoretical basis. 



With the appearance of detailed abu ndance Doppler im- 
ages (DI) of the stella r surfaces (e.g., [Luftinger e t alJ 12003b 
Kochukho v et al.l l2004t ILehmann et al.l 120071) and advanced 
computers it became possible to numerically explore the role of 
uneven distributi ons of chemical elem ents as a source of the light 
variability. First, Krticka ^t al.l (120071) succeeded in reproducing 
the light variability of a hot (r eff = 22000 K) B p star HP 37776 
based on the DI maps of He and S i derived by Khokhlova et al.l 
(I2000I) . Later, iKrticka et ail d2009h successfully fitted the light 
curves of a cooler (r e ff = 14500 K) Bp star HR 7224 using Fe 
and Si maps presented by Lehma nn et alJ (|2007). The main re- 
sult of both studies is a very good, parameter- free theoretical ex- 
planation of both the shapes and amplitudes of the observed light 
curves. The authors thus proved, on the one side, that an inho- 
mogeneous surface distribution of chemical elements accounts 
for most of the light variability in these stars, and, on the other 
side, they independently confirmed the results of DI mapping. In 
spite of significant progress in this direction, the role of abun- 
dance spots (and individual elements) in CP stars cooler than 
B-type remains unknown. 
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In this paper, we continue to present results of the light-curve 
modelling of CP stars, concentrating on the brightest A-type CP 
star s UMa (HD 1 12185, HR 4905). This star has been studied 
extensively in the past and many good photometric observations 
are available, including the data obtained by space photometry 
experiments. The surface distributions of seven elements (Ca, 
Cr, Fe, Mg, Mn, Sr, T i) were derived via the DI technique by 
iLiiftinger et alj d2003l) . providing the necessary basis for theo- 
retical modelling of the photometric variability. 



2. Observations 

The CP star s UMa is a well-known spectral, magnetic, and pho- 
tometric variable star. The Ca n K line and profuse lin es of Cr H 
vary opposi t ely in strength with a period of 5^0877 dGufhnick 
193 1[ 1 1934 IStmve&Hiltnelll943l: ISwenssonl \l944 iDeutsch 
1947). Photometric variations are observed to h a ve the same pe- 
riod, but to have double wave s dProvin|[l953l : iMusieloket ail 
1980; iPvper & Adelmanl 1 19851 etc.) with the maximum light 
occurring at Can minimum chosen to be the phase 0.00. Far 
ultraviolet OAO-2 spectrometer observations of s UMa indi- 
cate str ong lig ht variations in the anti-p hase with optical ones 
dMolnaril 1975r . iMallama & Molndfl977l) . 

Spectropolarimetry of s UMa has detected its very 
weak variable magnetic field whose maximum coin- 
cides with the optical light maximum dBorra & Lanstreetl 
1980: iBohlender & Landstreetl Il990t iDonati & Semell 11990: 



Wade et al.l !2000). 



Until now, almost all pa pers dealing with variations in the 
light of s UMa have used the Guthnick ( 1931b ephemeris 



HJD(Ca II min. int.) = 2426437.01 + 5*082,7 x E, 



(1) 



where E is the epoch. 

In this paper, we use light curves observed in different photo- 
metric systems. Stromgre m uvby measurements were taken from 
IPvper & Adelmanl d 19851) . and extended by 10-color medium- 
band Shem akha photometry wit h filters similar to standard 
uvby filters dMusielok et al.lfl980l) . B and V ob servatio n s clos e 
to the Johnson standard were performed by Provl^ dl953h . 
Unfortunately, no reliable data has been obtained in the broad- 
band U filter, the available UBV observations of ISrivastav3 
(1989) been ususable due to their excessive scatter. The same 
concerns available data in the Bj and Vt passbands from Tycho 
catalo g, while Hipparcos Hp photometric measurements (ESA 
119971) are well calibrated and were used here in the light curve 
modelling. 

We also used ph o tomet ric data derived bv lMolnarl d 19751) and 
Mall am a & Molnarl d 19771) using spectrometry of e UMa per- 
formed by t he ultraviolet satel lite OAO-2 (more about the project 
is given in ICode et alj[l9 70). We also attempted to use spec- 
troscopic observations obtained by the International Ultraviolet 
Explorer (IUE) and Copernicus space missions available via the 
Multimission Archive at STSc fl Finally, we used a dataset of 
photometric observations carried out by the Wide Field Infrared 
Explorer (WIRE) obtained during June-July 2000, although we 
had to rejected several par ts strongly affected by instrumental 
effects and instabilities (see Retter et al. 2004, for details). 



3. Methods 

3.1. Model atmospheres 

To perform the model atmosphere c alculations, we use d the 
most recent version of the LLmodels ( Shul yak et al.l l2004) stel- 
lar model atmosphere code. For all calculations, local thermody- 
namical equilibrium (LTE) an d plane-parallel geometry were as- 
sume d. The VALD database (Pis kunov et aljfl995l : iKupka et all 
1999) was used as a main source of the atomic line data for com- 
putation of the line absorption coefficient. The VALD compila- 
tion contains information about 66x 10 6 atomic transitions. Most 
of them come from the latest theoretical calculations performed 
by R. Kuruc20. 

The global magnetic field of the star is ver y weak (of polar 
intensity on the order of a few hund red Gauss (IDonati & S emel 
1990; Boh lender & Landstreetl 1990D ). which allows us to ignore 
its possible effect on the variability of the outgoing flux. 

3.2. Calculation of the light curve 

In previous attempt s to simulate t he light curves of H D 37776 
dKrticka et alJl2007h and HR 7224 dKrtickaet al.ll2009l) . the fol- 
lowing approach has been taken: 

1 . The construction of a grid of model atmospheres for a num- 
ber of abundances combinations that cover the parameter 
space provided by DI maps. 

2. Interpolation of specific intensities (or fluxes) from the grid 
onto a combination of abundances for all meshes visible at 
a given phase of rotation. The resulting flux is then obtained 
by surface integration of specific intensities (or fluxes with 
assumed limb-darkening law). 

This procedure is applied to every rotational phase and wave- 
length interval occupied by a given photometric filter. Although 
simple, this approach becomes challenging when modelling the 
light curves of s UMa. Indeed, in both cases of HD 37776 and 
HR 7224, only two elements were mapped (He and Si for the 
former; Fe and Si for the latter). To produce a smooth light 
curve, it was enough to compute a grid containing tens of models 
with different sets of abundances. Obviously, with the increas- 
ing number of mapped elements, the number of models needed 
to perform abundance interpolation can become comparable or 
even larger than the number of surface elements originally used 
in DI calculations. 

Seven mapped elements of s UMa would require time- 
consuming computations. We instead use a slightly different ap- 
proach. To reduce computational expenses, we decrease the res- 
olution of DI maps by interpolating them from the original grid 
of longitudes and latitudes onto a more sparse one. Then, for ev- 
ery pair of new latitude and longitude, we compute a model at- 
mosphere with individual abundances that represent a given sur- 
face mesh. Next, for every model (i.e. surface mesh), we com- 
pute the light curves in a particular photo metric filter wit h the 
help of modified computer codes taken from lKuruczl d!993[) and 
extend ed by Hipparcos photometry with passbands from Bessell 
(2000). The expected magnitude at a given phase c p h aS e is ob- 
tained by surface integration of individual fluxes F, from all vis- 
ible meshes 



Cphase = -2.5 log 



h I 



Fi cos 8iu(0i)dS 



visible surface 



(2) 



1 http://archive.stsci.edu/ 



http : //kurucz . harvard . edu 
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Table 1. Stellar parameters of s UMa. 



Effective temperature T e g 


9000 K 


Surface gravity log g (cgs) 


3.6 


Inclination angle i 


45° 


Abundance ranges 


from DI maps 


Ca 


[-6.149,-4.628] 


Cr 


[-8.270,-1.559] 


Fe 


[-6.529,-2.582] 


Mg 


[-5.686,-3.285] 


Mn 


[-6.281,-4.890] 


Sr 


[-11.515,-3.967] 


Ti 


[-8.478,-6.561] 



where 6 is the angle between the normal to the surface element 
and the line of sight, u(9) describes the limb darkening for a 
given filter, which we assume to be a quadratic law, i.e., 

u(6) = 1 -a(l -cos 8)-b{\ -cos 6) 2 , (3) 

and U is given by 

U= J u(Gi) cos QidS. (4) 

visible surface 

We adopted the same limb darkening coefficients regardless of 
the true surface elements. In particular, the quadrat ic coef f icients 
a and b in each color were taken from the tables of lClaretl (2000) 
for the solar abundance model. For the Hipparcos Hp filter, these 
coefficients were taken to be the average between the Johnson B 
and V filters. Since there is no definite filter curve for the WIRE 
star tracker, limb darkening coefficients for Johnson V filter were 
used instead (H. Bruntt, private communication). Because of 
this uncertainty, the WIRE photometry is presented in this pa- 
per mainly for illustrative purposes. 

The original resolution of DI maps of s UMa is 33 latitude 
and 68 longitude points equidistantly spaced between [0, 180] 
and [0, 360] spherical angles. To reduce the number of calcu- 
lations, we decreased the resolution of maps to 15 and 30 lat- 
itude and longitude points, respectively. As proven by numeri- 
cal tests, this decrease in resolution provides abundance maps 
smooth enough for an accurate representation of the light curves: 
we found no visible differences in light curves when computing 
all 2244 models at the original resolution, and 450 models from 
the reduced grid (see Sect, discussion for more details). Thus, in 
the remaining light curve modelling presented in this work we 
always used a lower number of 450 surface elements than the 
2244 original. 

Table [TJ summarizes the basic stellar parameter s adopted in 
the present study taken from Luftinge r et alJ (120031) . The abun- 
dances are given relative to the total number ot atoms, i.e. e e i = 

fog(AWNtotal)- 

4. Results 

4.1. Improved ephemeris 

Observational data used in the following analyzes of rotation- 
ally modulated variations of e UMa were obtained in the past 70 
years during which the star has revolved more than five thousand 
times. Although we do not know the true uncertainty in the de- 
terminati on of the canoni cal rotational period P = 5^0887 pub- 
lished bv lGuthnickl dl93ll) it can hardly be better than +0?0001 ! 



E 




-0.02 1 1 1 1 1 1 1 1 — 1 

1950 1960 1970 1980 1990 2000 2010 

Fig. 1. The dependence of the primary m aximum phas es derived 
from Eq.|7]on time. V - BV photometry (iProvinl 1953), full cir- 
cle - UV photometry from OAO-2 s atellite dMolnarfl 19751). O - 
10-color medium-band phot ometry dMusielok & Madeil ll988). 
□ - Hipparcos photometry JeSAI 1 19971) . A ubvv photometry 
of iPvper & AdelmarJ (119851) . the oblique line represents the 
Guthnick's ephemeris from Eq. dTJ, and the star on it the phase 
position to which are referred the spectroscopic maps we used. 



Table 2. Maximum phases according to individual sources of 
photometric data. A<p is the difference between new phases and 
phases evaluated according former ephemeris, see Eq. dHJ. 



Year 


max. phase 


A(p 


N 


source 


1952.3 


-0.004(10) 


0.006 


36 


Provin(1953) 


1971.0 


-0.000(3) 


0.024 


38 


Molnar (1975) 


1975.7 


-0.001(5) 


0.028 


179 


Musielok et al. (1980) 


1982.1 


0.012(7) 


0.035 


39 


PvDer & Adelman (1985) 


1991.4 


-0.006(6) 


0.044 


175 


ESA(1997) 


1995.2 




0.047 




Luftinger et al. (2003) 



This uncertainty manifests itself in the uncertainty of 0U in the 
rotational phase determination, which is for reliable studies un- 
acceptable. 

This compelled us to improve the s UMa rotational pe- 
riod. For this purpose, we used all available reliable photo- 
metric data, especially the 36 BV measurements of IProvinl 
(119531) . 38 synthetic magnit u des d erived from spectrometric 
data in 191 nm by iMolnarl (1 19731) . 179 individual observa- 
tions in the ten - color Shemakha medium-band system taken by 
iMusielok et al.l dl980l). 80 S t romgr een uvby measurements ob- 
tained bv lPvper & Adelmanl (Il985l) . andthe 314 B T V T , and 122 
Hp measu rements derived from data of the Hipparcos satellite 
dESAlll997l) . Thus, we had 769 individual measurements cover- 
ing the period of 1952-1993 (see Table |2|; unfortunately mea- 
surements from the present time were absent. Measurements 
taken in various narrow- and broad-band filters were divided into 
10 groups according to their effective wavelengths - see Table[3] 

We also used a very extensive data-set of WIRE photometry 
dRetter et alJ [2004). which was unfortunately unreliable in sev- 
eral ways: 

1 . We have been unable to interconnect the timing of its mea- 
surements with standard HJD timing. 

2. Because of their large scatter and instabilities, we omitted 
the first part of measurements and then one day of measure- 
ments. 
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rotational phase 

Fig. 2. Observed light curves constructed from original data 
in several passdands, whose effective wavelengths in nm are 
given on the right. For clarity, individual LCs are shifted in 
their ma gnitudes. Full circles - UV photometry from OAO-2 
satelli te (Molnai|| 19751) . O - 10-color mediu m-band photometry 
(Musi elok et al.lfl980l). A ubvv photometry of iPvper & Ade lman 
il985l) . V - fly dProvinll953l) . □ - Hipparcos photometry (IESaI 
1997), an d small full circl es relate to broad-band photometry of 
WIRE dRetteret all 120041) - each point is the median of 1600 
individual time-consecutive measurements. 



3. We found that consecutive measurements are not indepen- 
dent. We therefore aggregated WIRE measurements into 
groups of about 1600 members. Altogether we used 195 nor- 
mal points with standard uncertainty of 0.55 mmag. 

4. We identified significant trends in the observations, but 
which could be represented well by cubic polynomials and 
removed. 

The double-waved light curves are more or less similar (see 
Fig. |2]i, diff ering only in their effect ive amplitudes A c (for a def- 
inition, see Mikulasek et al. 2007a), where the subscript c de- 
notes the photometric passband. Light-curve magnitudes can 
then be expressed as 



m c j(f) — m c 



+ ±A C F(#), 



(5) 



where m C j(i) is the magnitude in color c observed by the y'-th 
observer, m C j is the mean magnitude, which can be variable over 
the long term. This was the case for WIRE measurements, where 
we had to assume a cubic trend. 



Table 3. Dependence of the effective amplitude of light varia- 
tions on wavelength. N denotes the number of measurements. 



^cff [nm] 


A eft [mag] 


N 


effective filter 


191 


0.1232(25) 


38 


Molnar 


350 


-0.0336(22) 


38 


u, U w 


360 


-0.0325(33) 


18 


Pw 


410 


-0.0147(27) 


39 


v, X w 


440 


-0.0190(22) 


175 




465 


-0.0218(18) 


38 


b, Y w 


510 


-0.0220(15) 


140 


H p , Zio 


550 


-0.0209(01) 


407 


V, y, V T , Vio, WIRE 


635 


-0.0230(33) 


37 


HRio, S io 


740 


-0.0205(28) 


34 


MR 10 , DRia 



-0.05 



CO 
CD 



CD 



CD 
> 



o 

CD 



0.05 



0.1 



1 0.15 



0.2 



+ 
+ 



2.2 



2.4 

l0 9irAfl 



2.6 
[nml 



2.8 



Fig. 3. Effective amplitudes of observed light curves taken in 
various effective wavelengths. Squares correspond to the ampli- 
tudes of LCs taken in i ndividual photometric filter s, full circles 
are values adopted from Mall ama & Molnarl (Il977l) . and crosses 
are effective amplitudes of LCs modelled in this work. 

Function F(§) is the simplest normalised periodic function 
that represents the observed photometric variations of s UMa in 
detail. The phase of maximum brightness is defined to be 0.00, 
and the effective amplitude is defined to be 1 .0. The function, be- 
ing the sum of three terms, is described by two dimensionless pa- 
rameters Pi and where Pi quantifies the difference in heights 
of the primary and the secondary maxima and Pi expresses any 
asymmetry in the light curve 



F(§, Pi , P 2 ) = yjl-P^-Pj COS(2 7T ff) + p X COS(4 7T ■&) + 
02 



-2= sin(2 n ■&) - 4= sin(4 n ff) 



(6) 



where # is the phase function. The O-C diagram (see Fig. Q] and 
Tabled indicates that the function is linear; thus, we assumed it 
to have the form 



# = (? - M ) IP, HJD(max T) = M + PxE, 



(7) 



where P is the period of the linear fit and Mo is the HJD time of 
the primary maximum nearest the weighted centre of all observa- 
tions except WIRE ones. The time of basic primary maximum in 
WIRE timing is then Mow- All 39 model parameters were com- 
puted simultaneously by a weighted non-linear LSM regression 
applied to the complete observational material. 
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We found these parameters to be P = 5^088631(18), M = 
2442 150.778(11), and Mow = 196 d 9647(31). The effective am- 
plitudes A c are given in Table[3j P\ = 0.557(5), and Pi = 
-0.056(7). The adequacy of the adopted linear model for the 
phase function given in Eq. (0 can be tested by studying the 
changes in the mutual phase LC maxima with time of observed 
light curves (see Fig. [I). 

Rotational phases evaluated according to Guthnick's 
ephemeris Eq. (Q}, ycuth can be transformed to new phases ac- 
cording to Eq. (O, y>new by means of the relation 



<£new = <£Guth + 2.648 x lO" 6 (f - 2 43 1 960) 
= VGuth + 9.670 x 10~ 4 (T - 1946.4), 



(8) 



where t denotes the time in JD, T the time in years, and their 
respectivefractions. 

Figure[3] represents the dependence of effective amplitudes 
of observed light curves on their effective wavelengths. Values 
summarized in Table[3] were calculated usin g UV effective am- 
plitudes derived by Mallama & M olnari d 19771) from spectromet- 
ric observations of the OAO-2 satellite. This figure clearly illus- 
trate the dependence of th e effective amplitude on wavelength. 
Mallama & M olnari (1 1 9771) predicted that the "zero point" should 
occur at a wavelength of 345 nm. However, the "zero point" (if 
any) appears to occur at a slightly shorter wavelength. 

Applying weighted LSM analysis to 28 measuremen t s of th e 
effective magnetic field p ublish e d inlBorra & Lanstreel dl 980). 
Bohlender & Landstreell d 19901) . iDonati & Semell dl990h .~~and 



Wade et al. (2000), we concluded that the observed sinusoidal 



variations can be satisfactorily reproduced by the simple model 
of centered magnetic dipole inclined by the angle f3 to the ro- 
tational axis (see Fig.|4j». The north magnetic pole then passes 
the central meridian at the phase <p = -0.002(21), and the am- 
plitude of the magnetic field variations is 142(17)G, their mean 
value being 24(6) G. We found the relation between angles p 
and the inclination of the rotational axis with respect to the ob- 
server ;', to be tan(/?) = 3.0(9) tan(/). Assuming (together with 
iLtiftinger et af]l2003h that i = 45°, we found that p = 72°(6°), 
thus the north magnetic pole is located close to the center of the 
leading photometric spot on s UMa, which is unlikely to be a 
coincidence. 

The magnetic field of s UMa is the weakest measured mag- 
netic field in magnetic CP stars which would be hardly de- 
tectable in any more rapidly rotating and fainter star. Its max- 
imum (polar) surface magnetic field is weaker than 400 Gauss. 
Thus we could neglect its influence in our computation. 

4.2. Abundance anomalies and emergent flux 

Before presenting the results of the light curve modelling, we at- 
tempt to verify model predictions about the impact of individual 
elements on overall energy distribution. Figure illustrates the 
synthetic energy distributions for models computed with maxi- 
mum, minimum, and mean abundances of every element from 
DI maps. This exercise illustrates the maximum possible energy 
redistribution effect caused by enhanced chemistry inside a spot. 
From Fig. [5] one can note two general features. First, there are 
only two elements that have (on average) the strongest effect on 
radiative energy balance: Fe and Cr. This is mainly because they 
have the largest number of strong lines in the UV and thus ef- 
fectively block radiation at those frequencies. Absorbed energy 
is then redistributed at optical wavelengths. This is illustrated 
by the positive delta's of uvby passbands shown on every plot. 
For Cr, the difference in the y-band between models computed 



150 - 



100 



50 



-50 



-100- 




-0.2 



0.2 0.4 0.6 0.8 
rotational phase 



1.2 



Fig. 4. Effective magnetic field and its uncertainty in gauss 
versus rotational phase. Note that the field is extraordinarily 
weak and the phase of its maximum coincides exactly with 
the maximum of brightness. Measurements taken form: o - 
iBorra & Lanstreetl d 19801). O - iBohlender & Landstreel d!990l) . 
A - IDonati & Semell d 1990) . □ - IWade et al l (l2000h . Areas of in- 
dividual markers correspond to the weight of the particular mea- 
surement, which is inversely proportional to the squares of their 
uncertainties. 



with mean and maximum Cr abundance is almost zero due to the 
strong Cr absorption features around 5500A that keep fluxes al- 
most unchanged for the highest abundance of [Cr] = -1.56 dex. 
This never happens for Fe, which has a smooth flux excess in 
uvby bands. Both elements contribute significantly to the 5 200 A 
flux d epression frequently found in CP stars (see iKupka et al.l 
|2003L and references therein). 

Models with enhanced Mg, Mn, and Sr cause only marginal 
or small changes in photometric filters. Interestingly and in con- 
trast to any other element, overabundant Sr leads to a noticeable 
flux deficiency in the v parameter because of some very strong 
Srn /14078, 4162, 42 16A lines, while uby passbands exhibit a far 
lower sensitivity to abundance changes (the same applies to the v 
filter in the Ca enhanced model because of the Ca n H & K lines, 
but with much smaller amplitude). This effect is even compara- 
ble to that of enhanced Fe. However, the mean abundance of Sr 
is only -9.6 and, on average, it is difficult to expect any notice- 
able impact of it on the resulting light curve. This is investigated 
in more detail later in this work. Finally, the contributions of 
Ca and Ti to the magnitudes of Stromgren photometry are very 
small. 

Although the relative roles of different elements can be esti- 
mated by the aforementioned simulations, to obtain quantitative 
results it is still necessary to carry out accurate modelling of light 
curves by taking into account the individual contributions of all 
spotted elements. 

4.3. Predicted light variation 

Using DI maps of individual elements, we verified the relative 
contribution of each of them to the total light variation. The 
result is presented in Fig. [6] We note that WIRE data were 
smoothed over 200 points to provide more or less a reasonable 
view, yet some instrumental effects are still visible. 

As expected from the previous paragraph, among the seven 
mapped elements there are only two that strongly contribute to 
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Fig. 6. Differential effect of inhomogeneous surface distributions of individual elements on light variations of s UMa. Fe - full 
thick, Cr - dashed , Ca - dash -dotted, Mg - dash-do t-dotted, Mn - lon g dashed, Sr - dott ed, Ti - thin dashed . Le ft panel: □ - 
iPvper & Adelmanl d 1985b . O - iMusielok etail {l980). Right panel: O - iProvinl (119531) . A - iRetteret all d2004 . □ - lESAl d!997l) . 
See online version for color figures. 



the amplitude of light curves in all presented photometric pa- 
rameters: Fe and Cr. Next follows Mn, whose contribution can 
only be recognized by a keen eye (mostly in b and V filters). 
The impact of the remaining elements is negligible. The cumu- 
lative impact of all seven elements on the light curves of s UMa 
is presented in Fig. [7] (thick solid line). A very good agreement 
in terms of both the shape and amplitude of the light variability 
is found for almost all photometric bands (see Fig. [3}. The only 
exceptions are passbands shortward of the Balmer jump, namely 
Stromgren's u, Shemakha U\q, and P\q for which the observed 
amplitude is approximately two times higher than the predicted 
one (see Fig . [3j ■ In principle, this may be a signature of som e 
other elements that were not mapped in Liifting er et alj (|2003), 
but still contribute to the light curve. Interestingly, a poorer fit to 
the u p assband than to others was also reported by Krticka et al. 
(2009) (see their Fig. 10). Nevertheless, the amplitudes in other 
passbands are well reproduced in both the present work and 
those cited above. Taking into account the large scatter in the 
observed points in general for the ubvy system, it is difficult to 
explain the observed discrepancy in M-band in terms of impact of 
some missing element, at least without additional DI mapping, if 
ever possible. Finally, Hipparcos photometry exhibits significant 
scatter with a peak amplitude that is a little larger than the one 
theoretically predicted. 

The forms of the simulated LCs agree with the observed ones 
very well, all of the modelled curves being systematically shifted 
with respect to the observed light curves by +01014. 

Figure |7] allows us to conclude that the inhomogeneous sur- 
face distribution of chemical elements, being derived spectro- 
scopically with the help of DI techniques, can quantitatively de- 



scribe the observed light variability detected in a number of in- 
dependent photometric systems. 

4.4. UV variations 

As can be seen from Fig. [5] the main source of light variability is 
the energy redistribution from UV to visual due to the enhanced 
abundances of certain elements. Thus, the fluxes in these two 
spectral regions should display an anti-phase behaviour. This 
has been experiment ally already reported for some CP stars (see 
ISokolovl2 000. 2006, for CU Vir and 56 Ari, respectively). Light 
curve modelling of HD 37776 and HR 7224 has also success- 
fully predicted this characteristic signature of stellar spots. 

For e UMa, we again had a problem with the lack of well- 
calibrated phase-resolved UV observations that could be used to 
study the effect of flux redistribution. The IUE archive mostly 
contains spectra obtained with a small aperture. This means that 
some of the stellar flux was missed during individual exposures, 
which is important to note in our study of the true variations in 
the continuum flux. Unfortunately, the only two spectra obtained 
with a large aperture (LWR06920RL and SWP07944RL) were 
obtained at the same rotational phase (<p = 0.166). 

Another possibility would be to use spectroscopic and pho- 
tometric data obtained by the Copernicus OAO-2 satellite. This 
mission was equipped with a number of photometers in the UV 
and visual . The descrip tion of the photometric passbands can be 
found in Mead^ (1 19991) . For instance, when analysing the spec- 
troscopic data of OAO-2 in the spectral region 1050 - 1830A, 
Molnai] ( 1 1975b detected significant flux variation caused by 
changes in line-blocking as a function of the rotational phase. 
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Fig. 7. Same as in Fig. [6] but taking into account surface distributions of all elements simultaneously. Theoretical predictions are 
shown for linear (dash-dot-dotted) and quadratic limb darkening laws based on the assumptions of [M/H]=+0.3 (dashed), [M/H]=0.0 
(full thick), and [M/H] =-0.3 (dash-dotted). 



An anti -phase variation in particular was reported for both spec- 
troscopic (1590 - 1690A) and photometric (OA O-2 /tl9 1 0A fil - 
ter) observations relative to the optical data of iProviiil (Il953l) . 
Unfortunately, the available online co-added spectroscopic data 
from OAO-2 spectrometers contains spectra obtained at differ- 
ent wavelength and, even worse, at different phases. We were 
thus unable to choose a unique spectral region that had been ob- 
served at a sufficiently large number of rotational phases (the 
case would be a narrow region around 1200A, but only 6 phases 
are available and with large errors). Nevertheless, to confirm the 
anti-phase behavior of the UV flux relative to the visual, we per- 
formed an exercise similar to those presented in Fig. 5 of lMolnarl 
(1 19751) . which illustrates the light curve of the OAO-2 /11910A 
photometer. The predicted light curve was computed based on 
model fluxes in the region 1770A < A < 2050A, whic h cor- 
respo nds to the FWHM of the OAO-2 S3F1 filter (see |Meade 
119991) . The theoretic al magn i tudes are shown in Fig. [8] The V- 
filter observations bv lProvirj (1 19531) are also plotted. The oppo- 
site variations in the optical and UV fluxes with respect to each 
other can be clearly seen. The double wave behavior of the sim- 
ulated O AO-2 curve is in excellent agreement with those pre- 
sented in Mornarl (ll975l) . yet the theoretical amplitude is smaller 
by a factor of two. We emphasize that there is a systematical de- 
viation between the predicted and observed magnitudes of the 
OAO-2 photometry, as illustrated in Fig. [3] Theoretical compu- 
tations display a large scatter blueward of log A = 2.4. At the 
same time, some predicted points do follow the observed linear 
trend of the dependence of amplitude on wavelength in UV re- 
gion. This is the case for filters at log A = 2.23, 2.4, 2.5. We note, 
however, that we did not attempt a detailed quantitative com- 
parison between the observed and modelled light curves or their 



amplitudes (also because the original OAO-2 filter data, as well 
as the filter values, are not available online), but only confirmed 
the anti-phase variations in UV and optical fluxes from the mod- 
elling point of view. 



4.5. Variations in the peculiar a-system 



The photometric parameter a is often used as a measure o f the 
peculiarity in CP stars. It was introduced by Mai tzenl (1 19761) and 
corresponds to the well-known flux depression at 5200A fre- 
quently se en in spectra of CP stars. Using mod el atmosphere cal- 
culatio ns, iKochukhov et alj ((2005) noted and iKhan & Shu lyak 
d2007l) later confirmed (on the basis of more detailed computa- 
tions) that Fe is the main contributor to the 5200A depression in 
the temperature range of CP stars. In addition, at low tempera- 
tures Si and Cr are also important. For HR 7224, Krticka et al. 
(2009) also identified the dominant role of Fe relative to Si. We 
carried out the same investigation and illustrate in Fig.|9]the im- 
pact of different elements on the variation in the a-index (note 
that, since the a-index is positive, the negative delta value in the 
plot correspond to the higher a at corresponding phases). We 
thus confirm the major role of Fe, although, the contribution of 
Cr is very important too. Moreover, at phase (p = 0.7 the con- 
tribution of Cr is comparable to that of Fe. A weak contribution 
of Ti can also be seen, but this fades relative to Fe and Cr. In 
general, the variation in the a-index appears to be small, on the 
order of a few mmag. A similar weak variat ion was reported for 
the hotter star HR 7224 dKrticka et al.ll2009l) . 
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clear visualization. 
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Fig. 5. Effect of individual elements on the synthetic energy 
distribution. The difference in fluxes in Stromgren uvby pass- 
bands between models computed with mean abundances and 
maximum abundance of the particular is presented on each plot. 
Fluxes were convolved with FWHM = 20A Gaussian for better 
view. The passbands (corresponding to the width of individual 
filter at half-maximum) of uvby filters are indicated by grayed 
feoxes. 



4.6. Variation in hydrogen lines 

We found a very small variation in the /? index o n the stan- 
dard system defined by [Crawford & Manderi ([1966) related to 
the Balmer Hf3 line, with a peak amplitude of « 1.8 mmag. 
This is an order of magnitud e weak er than the variations re- 
ported by iMusielo k & MadeJ (Il988l) . and may be the signa- 
ture of some additional mechanism affecting mainly profiles 
of hydrogen lines, as long as other photometric parameters 
are well fitted. For example, the rotational modulation of hy- 
drogen lines may be caused by a non-zero magnetic pres- 
sure in the stellar atmosphere, as reported by IShulvak et alJ 
(120101 120071) . An overview of oth er possible mechanisms can 
be found in IValvavin et alJ (120041) and we refer the interested 
reader to this work for more details . Add itional observations 
of the p index from IMusielok et al.l (Il980h have large errors 
and have significant scatter preventing the direct comparison 
with Musielok & Madei ( 1988). Last but not least, the data from 
IMusielok & Madei ( 1988) contains only three points at the max- 
imum of the [3 light curve and a definite outlier at if = 0.7, which 
should be interpreted with caution. We thus conclude that, if true, 
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the high amplitude of the observed variability in the B index can- 
not be attributed to abundance anomalies alone. 



5. Discussion 

s UMa is the third CP star for which a model-based light-curve 
analysis predicts excellent agreement between theory and ex- 
periment, based on the assumption of the inhomogeneous sur- 
face distribution of chemical elements derived by DI techniques. 
The same analysis o f two other an d hotter Bp stars H P 37776 
dKrticka et al.ll2007l) and HR 7224 dKrticka et al.ll2009l) . but us- 
ing slightly different methods, has also uniquely found that abun- 
dance spots may be a dominating mechanism of light variability 
in CP stars. 

The maximum differences between the predicted and ob- 
served light curves are seen in Stromgren u and Hipparcos Hp 
bands, as illustrated in Fig. [7] At the same time, other photo- 
metric parameters are fitted well enough to exclude any other 
element significantly influen c ing th e light curves. For instance, 
BV parameters from iProvinl dl953l) probably provide the clos- 
est fit, again, within the error bars of the observations (see three 
consecutive points at phase 1 .0, fo r instance). These er ror bars 
are larger for Stromgren passbands. Krticka et al. ( 2009) find the 
same larger discrepancy in u but for the hotter star HR 7224. This 
ensures that the impact of any missed element is less significant, 
yet not negligible. In principle, taking into acco unt the relative 
impor tance of different elements described in iKhan & Shulvakl 
(2007), Si may represent such a miss ed element, but it could not 
be mapped in [Luftinger et alj (120031) . The calibration of the u- 
band may also be problematic. Nevertheless, the remaining vby 
passbands can be fitted much better without the need to involve 
other elements in the modelling. 

As described above, to reduce computational time we de- 
creased the resolution of the original DI maps from 2244 sur- 
face elements to 450. The characteristic difference in computed 
magnitudes between these two sets of maps was found to be 
Ac « 6 x 10~ 4 mag, which is well below the detection limit. 

The choice of the limb darkening model used to reconstruct 
the light curves in different photometric bands can also influence 
the amplitude of variations. In our study, we used a quadratic 
law with coefficients for every photometric band (except those 
described above in the paper) taken from the solar abundance 
model. However, neither metallicity nor the choice of limb- 
darkening law seem to significantly affect the light curves, nor 
can they explain the discrepancy in u. This is illustrated in Fig. [7] 
where, as an example, we overplot the predictions of the lin- 
ear and quadratic limb darkening laws and the solar abundance 
model, as well as for the quadratic law with [M/H] = +0.3 and 
[M/H] = -0.3 models (note that the average metallicity derived 
from the DI maps is [M/H] = +0.02). 

6. Conclusions 

We have successfully simulated the light variability of Ap star 
s UMa caused by the inhomogeneous surface distribution of 
chemi cal elements detected by the DI maps of Luftinger et alJ 
J200l . A very close agreement in both the amplitude and shape 
of the light variability has been found for different and inde- 
pendent photometric systems. Our main conclusions are sum- 
marized below: 

- The presence of abundance spots on the stellar surface is the 
major contributor to the light variability. We did not intro- 



duce any free parameters to improve the agreement between 
theory and observations. 

- The light variability is due to the flux redistribution from the 
UV to visual region, which is initiated by the presence of 
abundance spots. 

- We support the conclusion of Khan & Shulyak (2007) that 
numerous lines of Fe and Cr are the main contributors to the 
well-known depression around 5200A. 

- The variation in the theoretical B index is very weak and 
approximately one order of magnitude smaller than the ob- 
served one. We thus conclude that inhomogeneous distribu- 
tion of chemical elements alone cannot explain the observed 
light curve in the B index. 

- The strong contribution of Fe and Cr to the flux variations 
allows us to conclude that, as for the hotter Bp stars that 
have been analysed so far (HD 37776, T s s = 22000 K and 
HR 7224, r eff = 14500 K) in cooler Ap stars it is also very 
likely that only a few elements significantly contribute to the 
light variation: Fe and Cr in the case of s UMa. This must 
of course be tested in stars with values of T e ff lower than 
r eff = 9000 K of s UMa. 
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